Time-dependent density functional theory of high excitations: To infinity, and beyond 



Mcta van Faassen^ and Kieron Burke & 
a Afdeling Theoretische Chemie, Scheikundig Laboratorium der Vrije Universiteit, 
De Boelelaan 1083, NL-1081 HV Amsterdam, The Netherlands 
b Departments of Chemistry and of Physics, University of California, Irvine, CA 92697, USA 

(Dated: February 11, 2009) 

We review the theoretical background for obtaining both quantum defects and scattering phase 
shifts from time-dependent density functional theory. The quantum defect on the negative energy 
side of the spectrum and the phase shift on the positive energy side merge continuously at E = 0, 
allowing both to be found by the same method. We illustrate with simple one-dimensional examples: 
the spherical well and the delta well potential. As an example of a real system, we study in 
detail elastic electron scattering from the He + ion. We show how the results are influenced by 
different approximations to the unknown components in (time-dependent) density functional theory: 
the ground state exchange-correlation potential and time-dependent kernel. We also revisit our 
previously obtained results for e-H scattering. Our results are remarkably accurate in many cases, 
but fail qualitatively in others. 
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I. INTRODUCTION 

Density functional theory (DFT) 0, H [| has become 
popular for calculating ground-state properties of large 
molecules; it replaces the interacting many-electron prob- 
lem with an effective single-particle problem that is com- 
putationally inexpensive to solve. Based on rigorous 
theorems [l[ll |i|, and a hierarchy of increasingly accurate 
approximations, DFT calculations are common in almost 
all branches of chemistry. 

In principle, all properties, including excitations, are 
functionals of the ground-state density, but accurate, 
general-purpose approximations to obtain these excita- 
tions have yet to be found. Instead, time-dependent 
density functional theory (TDDFT), based on another 
theorem [f|, provides a usefully accurate approach for 
many applications Q. In particular, low-lying single ex- 
citations are often well-approximated, and of great inter- 
est in many applications. Moreover, the computational 
technology for solving the necessary equations already 
existed, so TDDFT was rapidly built into existing codes 
and has become popular for calculating excitations. But 
TDDFT is mostly used to extract bound-bound transi- 
tions, and even then, only low-lying single excitations. In 
this perspective, we review our work extending TDDFT 
to higher excitations, including both Rydberg states and 
the continuum, i.e., n — > oo and beyond. 

We begin with a review of the theory of the quantum 
defect and scattering phase shift aimed primarily at our 
electronic-structure readers. Next, we give some back- 
ground on (TD)DFT and how to obtain the quantum 
defect and scattering phase shift from a (TD)DFT calcu- 
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lation. Wc prove that within TDDFT it is allowed to use 
the method of obtaining pseudo continuum states from a 
hard wall cavity, and thus obtain the phase shifts from a 
"bound" state calculation. We show results for the quan- 
tum defect of He and the e-He + scattering phase shift. 



II. THEORETICAL BACKGROUND 

The theoretical ingredients needed to understand this 
work originate from two distinct fields, scattering theory 
and density functional theory. While closely related in 
principle, these areas are almost mutually exclusive in 
practice. A reader from either field might skip the famil- 
iar material. 
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A. Scattering theory and Rydberg states 

The electronic structure of atoms is understood in a 
rather simple way if one of the electrons is in a highly 
excited state. In that case we consider the system as con- 
sisting of the excited electron in a local potential, orig- 
inating from the nucleus and all other (tightly) bound 
electrons. This local potential is a long-ranged Coulomb 
potential for neutral and positively charged atoms and a 
shorter ranged potential for negatively charged atoms. 

An atom with a highly excited electron in a long- 
ranged atomic potential is also known as a Rydberg atom. 
Such atoms are highly studied and an overview of their 
properties can for example be found in Refs. 0, H[ , and 
references therein. Since Rydberg atoms can be regarded 
as one electron in an effective local potential, they are 
similar to the hydrogen atom. The difference between 
Rydberg atoms and hydrogen atoms is expressed in terms 
of the quantum defect, which we introduce in the next 
section. 

The infinite Rydberg series merges into the continuum 
at the ionization limit of the system; it is this boundary 
between the finite states and the continuum that we ex- 
plore in detail in this paper. If the electron resides in a 
continuum state instead of a Rydberg state, it scatters 
from the resulting ionic core. We explain this further in 
section HTX2] 



FIG. 1: Accurate quantum defects for the He atom from 
Ref. [TO] (Exact) and experimental data from NIST [llj . 



I, and / is the ionization energy of the system studied. 
We use atomic units throughout. The quantum defect, 
/i;, is purely determined by that part of the potential 
that differs from —1/r. 

As is shown in Eq. (|2.I[) . the quantum defect behaves 
differently for different angular momentum states. For 
high angular momentum the quantum defect is gener- 
ally small, since the shielding effect of the core is mostly 
maintained. The existence of a quantum defect requires 
the atomic potential to decay like —1/r, producing a Ry- 
dberg series, but its value is determined by the interior 
region of the potential [{|. 

In Fig. [1] we show the exact quantum defects from high 
level wave function calculations [l(| , fitted to a smooth 
function as explained below, and we compare them with 
the experimental values [Tl| . We notice that both in the 
singlet and in the triplet case the quantum defect is a 
smooth function of energy, so we fit it to the following 
expansion, 

p 

^\E)=Y,^iE\ E = w-I. (2.2) 

i=0 



1. The quantum defect and its smoothness 

For any spherical one-electron potential that decays as 
— 1/r at large distances, the bound-state transitions form 
a Rydberg series. The high-energy levels of such a system 
are given by the Rydberg equation, which in atomic units 
reads, 
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(2.1) 



where i = (nl) and E l n is the excitation energy to a state 
of principal quantum number n and angular momentum 



In Fig.[fl the fits are given by n (1) {E) = 0.1395-0.0655£ 
for singlet and fi {1) (E) = 0.2965 - 0.0811S for triplet. In 
previous work [H, Gj, [HI we showed that in (TD)DFT 
this fit works well and wc only need two or three coef- 
ficients to describe all the data accurately. The reason 
for fitting the quantum defect is that it is hard to calcu- 
late the higher Rydberg states. If we have a way to fit 
the lower quantum defects accurately, we instantly have 
a way to obtain the values of all higher lying Rydberg 
states, up to the continuum. The coefficients for the he- 
lium atom can be found in the references above and are 
not repeated in this paper. 
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cross section is given by [lj 



(2.5) 



The scattering amplitude is written in terms of partial 
waves, 



f(e) = Y^f l p l (cosi 



(2.6) 



1=0 



where Pi (cos 6) is a Legcndrc polynomial and the partial 
amplitude can be expressed in terms of phase shifts: 
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FIG. 2: Helium atom negative energy singlet and triplet s- 
quantum defects from Ref. [Io( | and positive energy results 
from Ref. [l5|]. The phase shifts are divided by 7r so they 
match with the quantum defects at E = 0. The coefficients 
corresponding to the fits are reported in the text. 



2. Potential scattering and the phase shift 

The quantum defect, discussed in the previous section, 
is not only an important quantity for studying bound Ry- 
dberg states, but it is also relevant for scattering states. 
Seaton's theorem [l6| tells us that the quantum defect is 
a smooth function of energy as E — > 0, and it merges con- 
tinuously with the phase-shift (relative to pure Coulomb 
scattering) , 



5z(0) =tt W (0). 



(2.3) 



We illustrate this in Fig. [5] for the case of the helium 
atom. In this figure, and all that follow, we plot 7r/i; 
instead of the quantum defect itself. In this graph we 
also show a continuous fit of the points. Since the positive 
energy part is more curved than the QDs alone, we need 
more coefficients for an accurate fit. The fits are given 
by, 0.1399 - 0.0528S + 0.0358S 2 for the singlet case and 

0. 2967-0.0764£;+0.0167£; 2 for the triplet case. Together 
with these fits we also show the straight line fit of Fig. [T] 
It is clear that this fit is good for negative and very small 
energies, but diverges for larger energies. Eq. (|2.3p is 
only relevant if there is a quantum defect to consider, 

1. e. if we consider an atom with a long-ranged potential. 
Such a long-ranged atomic potential can in general be 
decomposed into a pure Coulomb (— 1/r) and a short- 
ranged part, 



v aT (r) 



v(r). 



(2.4) 



But to understand its scattering properties, we must first 
review elementary scattering from a spherical potential. 
If a potential is short-ranged, the differential scattering 



It follows from this equation that if we have the phase 
shift for all values of I we can calculate the differential 
cross section. In practice the few lowest I values will often 
be dominant and almost completely determine the cross 
section. 

For a long-ranged potential the equations are more 
complicated. In the case of atoms, we can separate the 
total potential into a Coulomb plus a short r ang e part. 
The differential cross section is then given by [17| 



/sr(M)| 5 



where fc(9) is the pure Coulomb amplitude and f av (9, </>) 
is the additional scattering amplitude as a result of the 
short-range part of the potential. This is necessary be- 
cause fc{@) diverges as 9 — > 0, making the cross-section 
itself divergent. Thus, scattering from a long-ranged 
potential is qualitatively different from that of a short- 
ranged potential. The additional scattering amplitude of 
Eq. (|2.8|) does not depend on the angle and can again 
be expanded as Eq. (|2.6[) . The partial wave amplitudes 
are given by, 
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sr,Z 
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(2.9) 



where the <t/ denote the Coulomb phase shifts and the 
8i the phase shifts due to the short-ranged part of the 
potential. If the short-range potential is small compared 
to the Coulomb potential, the additional amplitude can 
be calculated u sing , for example, a distorted wave Born 
approximation [171 ]. In the next section, we explain how 
we obtain accurate phase shifts in practice. 



3. Continuum states from spherical boxes 

A useful method we recently (re)discovered [HI, and 
which circumvents the problem of working with contin- 
uum states, is to put the system of scattering electron 
and target in a spherical box. In that way a finite num- 
ber of "continuum" states is obtained exactly, namely 
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FIG. 3: The lowest three positive energy orbitals for a spher- 
ical well within a box with Ri, = 10 (radius 1.2, well-depth 
3). The dashed curves indicate how the continuum functions 
continue. The curves are offset for clarity, each zero being 
marked by a dashed horizontal line. Color figure online. 



those that have nodes at the radius where we placed the 
wall, r = Rb. Assuming one has the liberty to choose any 
wall radius, all continuum states can be obtained in this 
way. Only for very low-energy continuum states does one 
need large boxes. 

Just as in the previous section, we distinguish two 
kinds of potentials, short-ranged potentials that effec- 
tively vanish after some finite radius R c , and long-ranged 
potentials that do not. The first kind we studied in 
Ref. where we investigated the scattering of an elec- 
tron by the hydrogen atom. The (Nt + l)-electron sys- 
tem, H - , is a negative ion with a short-ranged potential. 
Its phase shift is obtained from knowledge of the spheri- 
cal Besscl functions as a function of the energies and the 
wall radius, explicitly, 



tancL 



ji(k n R b ) 
ni(k n R b ) 



(2.10) 



where ji and ni are the spherical Bessel and Neumann 
functions. 

As an example, we show the lowest three of such states 
for a spherical well in Fig. [31 where the wall is placed at 
Rb = 10. The spherical well potential is given by 



V(r) = -DQ(a-r), 



(2-11) 



where D is the depth of the well, a the radius, and O the 
Heavisidc step function. The radius a of the spherical 
well is 1.2 bohr and the depth D = 3 Hartrcc. It is 
clear from this graph that the orbitals are the continuum 
functions, except that they do not continue beyond Rb = 
10 and that there is only a finite number of them. 

In this paper we also study electron scattering from 
He + . The (Nt + l)-electron system is the neutral helium 
atom, which has a long-ranged modified Coulomb poten- 
tial. We use the fact that this potential can be separated 
in a pure Coulomb potential (— 1/r) and a short-ranged 
potential, i.e. that vanished at some finite r. All we 
need to know to obtain the scattering phase shifts are 



FIG. 4: Same as Fig. [3] but with —1/r added to the potential. 
Color figure online. 



the Coulomb functions as a function of energy and wall 
radius, so Eq. ()2.10|) becomes 



tan 8 n 



Fi(r], k n R b ) 
Gi (rj, k n R b ) ' 



(2.12) 



where Fi and G; arc the regular and irregular Coulomb 
functions. The only requirement on the location of the 
wall is that it is outside the short-ranged part of the 
potential. We can place the wall relatively close to the 
atom, even though we are dealing with a long-ranged 
potential. 

As an example, we show the lowest three of such states 
for the same spherical well as above, plus a Coulomb 
potential, in Fig. [3J The total potential is given by, 



V(r) = -DQ(a-r) - 



1 



(2.13) 



Again the wall is placed at Rb = 10. It is again clear that 
the functions behave exactly like continuum functions. 



B. Ground-state density functional theory (DFT) 

We now review briefly both ground-state and time- 
dependent density functional theory We emphasize 
throughout the distinction between exact DFT and the 
approximations invariably used in practice. 

Density functional theory is based on the theorem by 
Hohenberg and Kohn (HK) Q. This theorem states 
that, for a given non-degenerate ground-state density of 
Fermions, n(r), with a given interaction, the external po- 
tential that produces the density is unique (up to an addi- 
tive constant). Thus, if the density is known, the external 
potential is also known, completely defining the Hamilto- 
nian, i.e., we can obtain all properties of the system from 
the density alone. The HK theorem also holds for de- 
generate ground-states [1| and there exists an analogous 
theorem that can be applied to spin densities [2(| Hl| . 

The power of DFT is in the idea of Kohn and Sham 
(KS) to apply the HK theorem to a system of non- 
interacting particles. The HK theorem states that the 
KS potential, which reproduces the true density of the 
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FIG. 5: Top panel: extremely accurate radial density for the 
helium atom found via the QMC method [22| • Bottom panel: 
The external and KS potentials for the helium atom. 

system, is unique. KS map the fully interacting problem 
to a non-interacting one with the same density. The KS 
equations are a set of one-particle Schrodingcr equations, 

(-^V 2 + v s (rf) fc(r) = EiMr), (2-14) 

where the 4>i are the KS orbitals and £, the corresponding 
orbital energies. The effective KS potential is given by 

v s (r) = v ext (v) + [ dv'^\ + v xc (r). (2.15) 
J |r-r'| 

The unknown part in this equation, u X c(r), is the 
exchange-correlation (xc) potential, that needs to be ap- 
proximated in practice. When the exact xc-potential is 
known, the exact interacting density of the system can 
be obtained from the KS orbitals, 

N 

n(r) = ]>> 4 (r)| 2 , (2.16) 

4=1 

where the sum runs over all occupied orbitals. 

We illustrate the KS idea in Fig. [5] where we show both 
the exact KS potential of the helium atom [22| and the 
external potential from an accurate wave function calcu- 
lation. Umrigar and Gonze [HJ obtained the exact KS 
potential from an extremely accurate quantum Monte- 
Carlo calculation for the ground-state of He, calculating 
the density, and finding the KS potential by inverting the 
KS equation for the doubly-occupied Is orbital: 

w B (r) = e u - V 2 <j> ls {r)/<f> ls (r), (2.17) 

where 4>\ s (r) = y^2n(r). The KS potential is quite differ- 
ent from the external potential. This is because two non- 
interacting "electrons" in the KS potential reproduce the 



FIG. 6: KS potentials for the helium atom: LDA, exchange 
only (EXX), and exact. 

exact electron density of the interacting system, by defi- 
nition. For example, because electron-electron repulsion 
pushes density outward, the KS potential is shallower 
than the external potential. 

The energy of the highest occupied molecular orbital 
(HOMO) of the exact He KS-potential is -24.592eV. 
This is exactly minus the ionization energy of the he- 
lium atom; an illustration of the exactness of Koopman's 
theorem in DFT [23[. This exact relation is a remark- 
able fact, since the KS system is a fictitious system, but 
it follows from the exactness of n(r) as r — ► oo. In fact, 
!) xc (r) — > — 1/r as r — > oo, for all spherical systems whose 
density n(r) decays exponentially. Thus, since, 

Mr) = / dr 1 - — - — — -h. - r -> oo, (2.18) 
J |r — r'| r 

then v s (r) — ► —1/r if N — Z , as in Fig.O or more rapidly 
if N = Z+1. 

In practice, v s (r) is not known exactly. DFT calcu- 
lations use approximations to Ey^ c [n\, whose derivative 
yields an approximation to u xc (r). Many approximate 
xc-encrgy functionals have been developed since the start 
of DFT. The simplest of these, the local density approx- 
imation (LDA), gives remarkably good results 0. It is 
improved by including the gradient of the density, lead- 
ing to the generalized gradient approximation (GGA). 
Some popular GGAs are BLYP (B88 for exchange 
and LYP [25] for correlation) and PBE (26|. Hybrid 
potentials are popular among chemists because of their 
improved accuracy for main-group thermochemistry and 
transition-state barriers. These hybrid functionals mix in 
some fraction of exact exchange with a GGA. The most 
widely used hybrid, B3LYP, contains 3 experimentally 
fitted parameters [25|, [2l], Hfl, whereas PBE is derived 
from general principles of quantum mechanics [26j . 

A less well-known feature to users of ground-state DFT 
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is that although their favourite approximations often 
yield very good energies (and therefore structures, vibra- 
tions, thermochemistry, etc.) and rather good densities, 
they have badly behaved potentials, at least far from nu- 
clei. The LDA potential, for example, decays much too 
fast. We show this in Fig. [6] for the helium atom. Al- 
though the relative energies of low lying states might be 
accurate, the high lying Rydbcrg states, in which we have 
particular interest, are poorly reproduced, if bound at all. 
In the case of the H~ anion, the LDA even fails to bind 
this system at all. These errors are called self-interaction, 
because E x = —E H and E=0 for one electron, but 
not with these approximations. Over the past decade, 
the technology for treating orbital-dependent function- 
als, such as the Fock integral (so-called exact exchange), 
has been developed, and such functionals help cure this 
problem [2^]. This is called the optimized effective po- 
tential method (OEP)[l)|,[3l|,[32j]. The OEP method can 
handle any orbital-dependent functional including simple 
exchange, and then produces an asymptotically correct 
v s (r). 



C. Time-dependent density functional theory 

Next we give a concise overview of the basics of time- 
dependent density functional theory (TDDFT), focussed 
on obtaining excitation energies. For a more in-depth 
review, we refer to Ref. @, and references therein. 



1. General formalism 

Ground-state DFT does not produce the exact tran- 
sition energies. The reason is that the KS orbital en- 
ergies do not correspond to the true bound-state en- 
ergies. A popular way to obtain the true transition 
energies is via time-dependent density-functional the- 
ory (TDDFT). TDDFT is an in principle exact method 
based on the Runge-Gross theorem Q. This theo- 
rem states that, for time-dependent problems, the time- 
dependent potential, and consequently all other proper- 
ties, is a functional of the time-dependent density. Armed 
with this theorem, time-dependent Kohn-Sham (TDKS) 
equations are defined analogues to Eq. (|2.14|) but with 
id(f>i/dt on the right. These equations describe a sys- 
tem of non-interacting electrons that evolve in a time- 
dependent Kohn-Sham potential and produce the same 
time-dependent density as that of the interacting system 
of interest, thus replacing the interacting problem with a 
much simpler non-interacting one, as in the ground-state 
case. One therefore defines (and needs to approximate) a 
time-dependent xc-potcntial, ti xc [n](r, f), as that defined 
by Eq. (12~T5D . 



2. Linear response 

Most practical TDDFT calculations use linear re- 
sponse theory. Within linear response theory, the 
TDDFT equations simply reduce to an eigenvalue equa- 
tion [33^ 

fi^Ff 5 = cufF^ (2.19) 

where the uji are the excitation energies of the system at 
hand. The superscript ± means that we have one set of 
equations for singlet energies, and one for triplet ener- 
gies. The key quantity in this equation is obviously the 
matrix, diagonalization of which, when spin-resolved 
for a spin-restricted system, immediately yields the tran- 
sition frequencies cof. 

n S& = 5 ab 5 t] {e a -e t f (2.20) 
+ 2V^Ti{KJl jb ±Kjl jb ) v/i^I" 

where the e 2 : are occupied Kohn-Sham (KS) orbital en- 
ergies and e a are unoccupied KS orbital energies, the K 
matrix elements are discussed below, ] and j specify the 
spin orientation. The Kohn-Sham orbital energies corre- 
spond to a system of non-interacting electrons, the sum 
of which, in an exact DFT calculation, gives the true 
density of the system. The first term of the is the 

KS transition energy. The second term in both matrices 
corrects the KS result to give the true interacting result, 
the singlet and triplet transition energies respectively 

An approximation to solving the above eigenvalue 
equation is the small matrix approximation (SMA), 
where we o nly include diagonal elements of the fi^- 
matrix [34|, [35[ ■ This leads to a singlet triplet splitting 
but it does not allow different transitions to mix with 
each other. The SMA is useful for understanding results, 
but usually is not applied to obtain the transition ener- 
gies for practical systems. 

We illustrate how the transition energies evolve from 
KS to SMA to full TDDFT in Fig. □ where we consider 
the Is — > 2s and Is — ► 3s transitions in the helium atom. 
From a ground-state KS calculation we obtain just the 
KS orbital energy differences, which are somewhere in- 
between the exact values. The next step is to use the 
small matrix approximation (SMA), where we see that 
especially the singlet values still overestimate the exact 
values. The final step is to include the full matrix, al- 
lowing full mixing of all states. The full TDDFT results 
are close to the exact values, but do not match them, 
because we used an approximate kernel (see below). Of 
course, all this information, and much more, is beauti- 
fully condensed into the quantum-defect plots through- 
out this paper. For example, the exact results in Fig. [7] 
are represented by the leftmost points of Fig. [TJ 
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FIG. 7: The first two transitions for He. Exact KS orbital 
energy differences, SMA singlet-triplet excitation energies and 
full TDDFT (ALDA) excitation energies (both evaluated for 
the exact ground-state KS potential), and the exact values. 



3. XC Kernels 



The K matrix elements in Eq. (|2.20p are given by, 
<ka Q)(/w (r)0b T (r')</> JT (r') 



r - f 



^(r)</w(r)/ x 7(r, r', ^)^ T (r')^.(r') drdr(2.21) 



where we used er, r to indicate spins and i,j indicate oc- 
cupied orbitals and a, b virtual orbitals. The first term in 
this equation is the Hartree term, which is only depen- 
dent on the ground-state KS orbitals. From Eq. (|2.20p 
we see that this term cancels in the triplet case (minus 
combination). The second term contains the xc-kernel, 



5V XC< y(V,t) 

<W(r',t') 



(2.22) 



whose Fourier transform is simply /^ c T (r, r', lo). The xc- 
kernel is the only part in this expression that needs to 
be approximated, after the ground-state calculation has 
been performed. In this paper we focus on three differ- 
ent kinds of xc-kernels, adiabatic LDA (ALDA), EXX, 
and hybrid, and study how they perform for the bound 
and continuum states of the helium atom. All these ker- 
nels use the adiabatic approximation, i.e., ignoring the 
frequency dependence of the xc-kernel. 

The simplest xc-kernel is the adiabatic local density 
approximation (ALDA). The ALDA xc-kernel is given 
by 



/•ALDA 



T[n }(v,v',t,t') = 
<5 (3) (r - r') S(t - t') 



dn„dn a i 



, (2.23) 



n<,=n 0<T (,r) 



frequency-dependence and via a Kramers-Kronig rela- 
tion, this implies that it is purely real (3fj| . 

Even though finite systems like atoms and molecules 
are far from representing the homogeneous electron gas, 
the ALDA is a good approximation in practice. There 
are systems, though, where the ALDA fails completely. 
For example the polarizability of long chain molecules 
is strongly overestimated by the ALDA [37| ■ A solu- 
tion to this problem is to use orbital dependent func- 
tionals [38l [39j . or to solve the time-dependent cur- 
rent DFT equations with the current dependent Vignale- 
Kohn functional [4?j| . 

In this paper we study the helium atom, in the special 
case of a two-electron system we can make some other 
simple choices for the xc-kernel. For two-electrons we 
can treat exchange exactly, the kernel is simply minus 
half of the Hartree kernel, 



EXX 
xtrtr' 



Sera' 



r - r' 



(2.24) 



This kernel can be improved by combining with the corre- 
lation part of ALDA. This results in the so-called hybrid 
xc-kernel, 



j-hyb 
J XCcrcr' 



1 s„ 



2 r - r' 



1 f ALDA 
2 •/ caa' i 



(2.25) 



which improves on EXX for the helium atom [4l| . 



III. ONE-DIMENSIONAL ILLUSTRATIONS 

In this section we treat potentials of finite range in 
one dimension, i.e., v(x) = beyond some a. We re- 
strict ourselves to symmetric potentials (v(— x) = v(x)), 
where the analogy to 3d scattering theory is most plain. 
These examples are not meant as models of real atoms 
or molecules, but simply illustrate the general principles. 



A. Potential scattering: Nt = 

We divide the illustrations into two cases, to make as 
clear as possible the distinct physical principles involved. 
We first treat potential scattering, i.e., a single-particle 
scattering from a potential, and show how box calcula- 
tions yield the exact result. Thus Nt = 0, where Nt is 
the number of particles in the target. In the next section, 
we treat the more complex interacting case, analogous to 
electron scattering from H, where Nt = 1. 



1. Phase shifts in Id 



where e"" lf is the xc-density of the uniform electron For a symmetric Id potential, parity plays a role analo- 
gas. The time Fourier-transform of this kernel has no gous to the partial waves in 3d. Thus, there are only two 
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channels, symmetric and antisymmetric. Careful analy- 
sis allows definition of two phase shifts 



e 2iS± =t±r. 



(3.26) 



The + sign corresponds to an even continuum solution 
and the - sign to an odd solution. The amplitudes t and 
r are defined for scattering boundary conditions, i.e., the 
asymptotic wave functions are, 



(j)k{x) — ► e lkx + re lka 
-> te ikx 



x — ► — oo 

x — > oo. (3.27) 



The amplitudes r and t are related to the reflection and 
transmission coefficients, R and T, 

R=\r\ 2 , T=\t\ 2 . (3.28) 

A useful result from this analysis is 

i = e i05 + +<n cos A, (3.29) 

where A = 8 + — 8~ , showing that resonances occur when 
A /it is an integer. 

An example of a short-ranged central potential is the 
delta well potential, v(x) = —ZS(x). The free space ex- 
ample can be easily calculated, and the amplitudes t and 
r are given by 



t 



1 



1 + Z/ik' 

From this it follows that 

2*5+ = 1 - Z/ik 
1 + Z/ik' 

and we find for the phase shift 

r+ Z 

taii(V = — . 



r = t- 1. 



(3.30) 



(3.31) 



(3.32) 



The odd phase shift S =0. (This is special to the delta 
potential, but is not true in general.) 



2. Scattering solution in a box 

Now consider putting the potential inside an infinite 
well, with walls located at x = ±L/2. Since the poten- 
tial is placed in a box, the states are discrete and there is 
only a discrete set of allowed values for the eigen-energies 
E„, n = 1,2, . . .. By comparing wavefunctions beyond 
the range of the potential, one finds the phase shifts cor- 
responding to the even and odd solutions are given by 



4 2 4' 



(3.33) 




FIG. 8: Solid line: even continuum state for a delta function 
with strength Z = 2 at k = 4.2748, corresponding to the first 
state of positive energy when walls with L = 2 are imposed. 
Symbols: bound state for the same potential at the same 
energy in a box with walls at x = ±L/2 = ±1. 



The even solution to the Schrodinger equation in a box 
with boundaries at x = ±L/2 is 



<t>ti x ) = -4sinfc+ 
with normalization constant 



A= 12 f sinfc+L \ 
V L \ k+L J 

The discrete allowed fc+ values satisfy 

k+ T, 

k + =Ztan^± 



-1/2 



(3.34) 



(3.35) 



(3.36) 



(For the derivation, see for example Ref. [43j|.) Taking 
the tangent of Eqs. (|3.33[) . wc obtain for the phase shift, 



tan S 



+ _ 



cot ■ 



k„ L 



Z 



(3.37) 



Eqs. (|3~32|) and (|3~37)l are identical, but Eq. |3~37)) can 

only be evaluated at a finite number of energies fc+. In 
Fig. [8] we show both a continuum scattering state and a 
box state for the delta well potential at the same energy. 
From the figure it is clear that they are identical within 
the box, and therefore have the same phase shift. 
The odd solution is 



sin k„ x. 



(3.38) 



with kz. 



\f2E^ 



These are simply the particle-in-a-box wave functions; 
the delta potential has no effect on the odd states. The 
corresponding energies are the particle in a box energies, 

2n7r 

K = — 

where n = 1,2,3,.... This corresponds to a phase shift 
5~ = 0. 



(3.39) 
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3. DFT treatment 

It is now instructive to see how this one-electron exam- 
ple should be treated with (TD)DFT. The first step in 
any linear-response TDDFT calculation is to find the KS 
ground-state potential of the (Nt + l)-electron system. 
The ground-state KS potential is given by 

vks{x) = v ext (x) + vnxc(x)- (3.40) 

For one electron (N = 1), v x (r) = — v H (r) and v=0, and 
the KS potential is just the external potential. The KS 
orbitals for positive energies are equal to the scattering 
states we obtained above in the exact case (with or with- 
out the box). 

Next one usually calculates the TDDFT correction cor- 
responding to the excitation of the electron into one of 
the positive energy states. Again, for a one-electron sys- 
tem, v H xc(r, t) = 0, and thus /hxc(i", t',uj) = 0. From 
Eqs. (pT2T))) and ([2l>Tj) it follows that the true excita- 
tion energies are equal to the KS orbital energy differ- 
ences. The TDDFT scattering phase shift is therefore 
equal to the exact phase shift in this trivial case. This 
also tells us the condition needed for any TDDFT treat- 
ment to recover the exact results for this case: Both the 
ground-state potential and the TDDFT kernel must be 
self-interaction free for one electron. 




FIG. 9: Real and imaginary parts of the exact (solid lines) 
and KS (dashed lines) transmission amplitude for the delta 
potential well with contact interacting electrons. In this plot 
Z = 2 and A = 0.5. 



identical to Eq. (|3.30p . But the singlet has a correction 
due to the interaction: 



t ■ — r 

^sing — L 



1 - 



2iXk 
Z 2 + k 2 



(3.43) 



Our methods should be able to reproduce this leading 
correction. 



B. 



Interacting case: Nt = 1 



The case in the previous section was a trivial one for 
(TD)DFT, but checked that our methodology worked in 
that limit. Now we study a case that is not trivial. We 
consider two fermions, one of which resides in the bound 
state of the delta potential well and the other scatters 
from this target. The two electrons interact through a 
delta potential repulsion. We show in this section that, 
with exact functionals, the TDDFT results reproduce the 
exact results. 



Exact solution 



The Hamiltonian of this problem is given by [44 



H 



Z [5{xi) + 8(x 2 )\ + X8(x 1 - x 2 ) 



2 \dx\ 

(3.41) 

The two electrons interact via a delta-function repulsion 
potential, scaled by A. In Ref. [45[ the scattering ampli- 
tudes for both the exact and (TD)DFT case were derived 
to first-order in A. They found the transmission ampli- 
tude for the exact case [45|. For triplet scattering, the 
bound electron has no effect, and 



t — t" — 

''trip — <- 



Z + ik 



(3.42) 



2. Exact ground-state KS potential 

To reproduce this result with TDDFT, first the 
ground-state KS problem must be solved. The ground- 
state KS potential of the (7Vr+l)-electron system is given 
by 



v s (x) = -Z5(x) + \\<p s (x){< 



(3.44) 



to first-order in A. It is possible to solve the correspond 
ing self-consistent Schrodinger equation exactly, finding 

<j> s (x) = VZe- z ^ + ^= \ze~ 3Z ^ + e~ z ^(4Z\x\ - 3)' 



(3.45) 



To first order in A, the KS potential is therefore 

v s (x) = -ZS(x) + \Ze~ 2Z W. (3.46) 

Transmission through this potential can easily be cal- 
culated using the distorted-wave Born approximation, 
which is exact to leadingorder in A, producing the KS 
transmission amplitude 451 ] . 



f 1 



i A k 
Z 2 + k 2 



(3.47) 



The transmission through the KS potential is exactly 
midway between singlet and triplet scattering. In Fig. [51 
we plot both the exact and KS transmission functions for 
Z = 2 and A = 1/2 to first order in A. 
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3. TDDFT corrections 

The next step in a TDDFT calculation is to correct the 
KS transmission for both singlet and triplet scattering, 
using the Hartrcc-exchange (Hx) kernel as correlation ef- 
fects are higher order in A. Because we are treating both 
singlet and triplet scattering, we must use spin-TDDFT. 
This has no effect on the closed shell ground-state, but 
the exact spin-dependent kernel is, to first order in A, [34[ 



(J(J 
HX 



\5{x-x'){l-S aa ,) 



(3.48) 



where a indicates the spin. 

Using this kernel the following amplitudes are ob- 
tained, by applying a distorted-wave Born approximation 
directly to the Dyson-like equation relating the true and 
KS response functions [45| : 



t, - X- 



'(k-iZ) 2 (k 



iZ) 



^trip — t s + A 



(k-iZ) 2 (k + iZ) 



(3.49) 



(3.50) 



These equations are identical to the exact equations 13.431 
and 13121 This illustrates that TDDFT reproduces the 
exact answer. 



1 in a box 



We are finally ready to perform the most important 
test of our methodology. We use TDDFT with exact 
ground-state and time-dependent functionals and with 
box boundary conditions to reproduce the exact inter- 
acting scattering results above. 



1. KS ground-state in box 

To begin, we need first solve the ground-state KS prob- 
lem in the box. Now, the KS potential, Eq. (|3.46[) . is 
not zero beyond some radius R, but extends out to oo. 
But since it is exponentially decaying, we can choose 
L » 1/Z, and have only exponentially small corrections 
to the free results. Then the KS potential inside the box 
is simply given by Eq. (|3.46p . For the even states, we can 
use the eigenstates Eq. 13.341 to zero-order in A, and then 
apply first-order perturbation theory. The zeroth-order 
eigenstates simplify to, 



4>t ( x ) = \l 7 sinfc^ 



(3.51) 



i.e., the normalization constant of Eq. 13.351 simplifies, by 
using the allowed values of fc+ and the condition L » 



1/Z. Note there is no assumption concerning the relation 
between fc+ and L. 

We next treat the interaction potential as a perturba- 
tion. The energies of the scattering states in the box can 
be expressed as, 



Se n 



(3.52) 



where e„ are the energies of the delta potential well with- 
out interaction. The perturbed energy is given by 

fL/2 

Se n = (n\ 5v \n) = 2XZ / dxe- 2Z ^\<j) n {x)\ 2 . (3.53) 
Jo 

An elementary integration yields: 



8e T 



1 



A 



L l + Z 2 /k 2 



(3.54) 



for either odd or even states, using Eq. (|3.35|) and L >> 
1/Z. Now, since 5k = 5e/k, 



L5k„ 



L 



Sv r , 
kn 



A 



1 + Z 2 /k 2 • 



(3.55) 



Finally, inserting these into Eqs. (|3.33[) for the corrections 
to the phase shifts, and then calculating the correspond- 
ing transmission amplitude from Eq. (|3.29[) . we obtain 
for the KS transmission amplitude for two contact inter- 
acting electrons in a delta well potential, 



t s = fO e -iVMi+zVtf) =i° [ 1 



i\kr, 



(3.56) 



Z 2 + ki 

agreeing with the exact result, Eq. (|3.47j) . 

2. TDDFT corrections in box 



Last, we need to show that the TDDFT corrections 
for the box yield the exact results. We start from the 
^-matrix of Eqs. (|2T2U)l and (gZZTP . We first make the 
small matrix approximation (SMA) that assumes that 
only diagonal elements of this matrix contribute to the 
excitation energies. For the singlet excitations, we can 
then write the il-matrix in short-hand notation, 



u q + (q\fnx\q), 



(3.57) 



where q = (0, n) indicating that we only consider transi- 



tions from the ground-state. Writing , 
/ is the ionization potential, we find: 



2(<7|/ H xk), 



where 



<<?!/, 



x { q) = 2XZ I 
o 



L/2 



-2Z\x\ 



Q„ — /, where 



(3.58) 



\<l> n {x)\ 2 da, (3.59) 
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FIG. 10: Real and imaginary parts of the KS transmission 
amplitude and of the interacting singlet and triplet ampli- 
tudes for the delta potential well with contact interacting 
electrons. In this plot Z = 2 and A = 0.5. For the results 
from the box (symbols) we used L = 5. 
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FIG. 11: s-Wave phase shift for the He atom. Accurate data 
from Refs. 0, HE]- All data below the KS result corresponds 
to singlet values, all data above to triplet values. For E < 0, 
we plot 7r times the QD, so it merges smoothly with the phase 
shift. 



i.e., the same as 8e n from Eq. (|3.53p . We then immedi- 
ately obtain for the transmission coefficient, 

= '° " #^5) ■ < 3 ' 60 > 

The same as in the continuum case Eq. Q3.49p . In the 
triplet case / HX = and the equations we obtain are 
identical to the continuum case, i.e. itrip = t ■ 

Thus exact TDDFT in the box precisely reproduces 
both the TDDFT results without the box, and the exact 
results from wavefunction theory, as shown in Fig. 1101 



energies from the DFT calculation. The number of states 
included in the Sl-matrix are always chosen large enough 
for the transition energies to have reached convergence. 

Since we want to asses the performance of TDDFT ker- 
nels, we use a highly accurate ground-state xc-potential 
for the helium atom, obtained from very accurate wave 
functions, by Umrigar et al. j22j. We regard this as the 
essentially exact xc-potential. This potential is spline- 
fitted onto the grid points and imported into the code. 



IV. COMPUTATIONAL DETAILS 

Here we describe the details of the numerical calcula- 
tions for real atoms and ions. We obtain the positive KS 
orbital energies (necessary to evaluate the phase shift in 
Eq. (|2.12[1 ) from a well-established fully numerical spheri- 
cal DFT code (46j . This code includes many approximate 
xc-potcntials, including the standard LDA and GGAs, 
and also the optimized effective potential method (OEP) 
for exact exchange calculations. The program does not 
use a basis set like most quantum chemistry codes, but 
instead works with a radial logarithmic grid. Both the 
KS orbital energies and the potentials are optimized in a 
self-consistent way. The code is supplemented by the op- 
tion to insert a hard-wall at a distance Rb from the origin. 
We always use a large number of grid points, > 1500, to 
ensure convergence with the number of grid points. 

The TDDFT excitation energies are calculated with 
our own code that solves the radial TDDFT equations. 
Since we are dealing with small systems, we can exactly 
diagonalizc the fi-matrix of Eq. (|2.20p without regard to 
computational cost. The TDDFT code reads in the grid 
points, ground-state orbitals, and ground-state orbital 



V. RESULTS 

In this section we show our results for c-Hc + and e- 
H scattering. In case of He we also study the quan- 
tum defect. We study in detail the effect of different 
xc-potentials and xc-kernels. We also study the effect of 
the location of the potential wall. 



A. Usefulness of the KS potential 

In Fig. [IT] we show highly accurate s-QD values (mul- 
tiplied by 7r in this and all following figures), and s- 
wave phase shifts, S s obtained from the literature [lOL fl5l] - 
We note, for our DFT readers, how smooth these curves 
are, demonstrating the continuity of cigenstatcs through 
E = 0. All cigenencrgies of the infinite Rydbcrg series 
can therefore be captured by a single smooth curve, as 
was already indicated in Eq. (|2.2p . The fit coefficients 
corresponding to the continuous curves in this and all 
following figures are available as supplementary material. 
Comparing our KS results to the literature data, we see 
that the KS QD and phase shifts are in-between the ac- 
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FIG. 12: Error of calculated singlet s-wave phase shifts with 
respect to the reference data from Refs. [Tol. [l5j. For E < 0, 
we plot 7r times the QD error, so it merges smoothly with the 
phase shift error. Color figure online. 



curate wave function results. The KS results also mimic 
the correct shape of the accurate curves. We note for our 
non-DFT readers how much information is contained in 
the exact ground-state KS potential of the (N T + 1)- 
electron system (in this case He). All we need TDDFT 
to do is to shift the curve down in case of the singlet and 
upward in case of the triplet. Thus, an important mes- 
sage of this review for non-DFT readers is just how good 
a starting point the accurate KS potential is. 



B. TDDFT corrections 

To see how, and if, TDDFT corrects the KS results, 
we study the effect of the f xc kernel that appears in 
Eq. (|2.2ip by using a sequence of approximations that 
include more and more physics. While there is as yet no 
explicit reverse-calculation for the exact kernel, we know 
that, if we had that exact kernel, it would reproduce the 
reference curves of Fig. [TTJ 

In Figs. [T2] and Q1J] we show the difference of the 
TDDFT results with the exact results for the s-wave 
quantum defect and phase shift of singlet and triplet 
respectively. The difference is calculated by fitting the 
results to a polynomial of order 2 or higher, and sub- 
tracting the exact values from the TDDFT values for the 
quantum defect and phase shift. In case of the quantum 
defect we indicate with circles the locations of the finite 
Rydberg states. We note again that for fitting just the 
QD, i.e. obtaining the /i; coefficients, a polynomial of 
order 2 is enough, while we generally need a higher order 
polynomial to fit the QD and phase shifts together. We 
now analyze our results in much detail. 

The first step in our analysis is to set the xc-kernel to 



FIG. 13: Same as Fig. 1121 but for triplet. Color figure online. 



zero. In case of the triplet, Eq. (|2.2ip tells us that we 
are left with just the KS orbital energy differences, i.e. 
there is no correction. For the singlet excitations there 
is a correction. The approximation of setting the xc- 
kcrncl to zero is called the Hartrcc approximation and is 
also known as the random phase approximation (RPA) in 
physics and does not involve any TDDFT. Fig. [T2"l shows 
that the Coulomb kernel corrects the KS result in the 
right direction, but overcorrects by about as much as the 
original error. We observe that even without an xc-kcrncl 
this methodology splits the bare KS results into separate 
singlet and triplet curves. 

The next step is to include exchange, which usually 
dominates over correlation. Since we are considering a 
two-electron system, the exact exchange kernel (EXX) is 
also known; it is given by Eq. (|2.24p . We see in Figs [T2l 
and [T3l that the Hartrcc results are considerably corrected 
by inclusion of the exchange kernels. In the triplet case 
the EXX correction is in the right direction, but it over- 
corrects. 

Correlation is included by considering the full ALDA 
kernel. The ALDA kernel performs slightly worse than 
the EXX kernel in case of the singlet results. In case of 
the triplet results the ALDA performs better than EXX, 
especially at lower energies. The size of the error in the 
singlet and triplet full ALDA curves compared to the 
reference data is very similar. 

For completeness we also show ALDAx results, which 
includes only the exchange part of the ALDA kernel. In 
the singlet case, the full ALDA results merge with AL- 
DAx for higher energies, the full ALDA improving on 
ALDAx for lower energies. In the triplet case the AL- 
DAx results are better than full ALDA for all energies. 

The good results of the EXX kernel and the correla- 
tion correction of ALDA, suggests to combine the two. 
This leads to the hybrid kernel [H[ of Eq. ([2~25]) . In the 
singlet case the hybrid kernel clearly improves the EXX 
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and ALDA kernels over the entire energy range. One of 
the conclusions of Ref. [4l[ was that the hybrid kernel 
improves the (bound-bound) excitation energies of He in 
almost all cases. On the continuum side, we see that 
the hybrid kernel performs the best of all these kernels 
in case of the singlet, and for the triplet it improves on 
EXX, but only improves the ALDA for large energies. 
Another conclusion of Ref. [4l[ is that the hybrid kernel 
improves the bound-bound triplet s excitation energies 
of He. Our finding that the ALDA performs better in- 
stead is most probably because the TDDFT calculation 
of Ref. [4l| only includes the lowest 34 unoccupied states 
of s and p symmetry in the calculation of the fi-matrix, 
whereas we included over 350 states including the con- 
tinuum in all cases, notably improving the numerical ac- 
curacy of the TDDFT results. 



C. Approximate ground-state potentials 
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FIG. 14: The error in the s-wave phase shifts obtained from 
the EXX potential, relative to those from the exact KS po- 
tential. Color figure online. 



Until now we have always used the exact ground- 
state potential for our calculations. It is interesting to 
study the influence of the ground-state potential on the 
TDDFT results as, in practice, only approximate poten- 
tials are available. There is a wealth of such potentials 
and wc cannot cover all of them in this paper. 

For illustrative purposes we therefore limit ourselves to 
two potentials: the most simple potential of all, the local 
density approximation (LDA), and the orbital-dependent 
exact exchange potential implemented within the opti- 
mized effective potential scheme (OEP). We expect little 
qualitative change moving from LDA to GGA, or adding 
GGA correlation to EXX. We showed these potentials to- 
gether with the exact potential in Fig. [6] From the figure 
we see that the LDA is very different from the exact po- 
tential. The potential is too narrow in the interior region, 
resulting in a bound state of too high energy, and it is 
short-ranged, completely missing the — 1/r tail. But the 
shape of the LDA in the interior region is typically very 
good. It is almost equal to the exact potential, shifted 
upward by a constant, therefore occupied LDA orbitals 
yield a remarkably accurate density, except in the tail 
region. This feature can be used to extract a highly ac- 
curate QD from the LDA potential [47j . 

On the other hand, the EXX potential is almost indis- 
tinguishable from the exact potential; it is very similar 
in the interior region and has the correct asymptotic be- 
havior. The exact potential is slightly deeper than the 
OEP potential near the nucleus (see the inset in Fig.HJJ). 
We expect the results from the OEP potential to be very 
close to the exact results. 

We show the error in the phase shift obtained using the 
EXX ground-state potential in Fig. [H] (using the ALDA 
kernel) . The LDA results diverge completely and there is 
no quantum defect because of the missing Rydberg series, 
therefore we do not show these results in the figure. The 
EXX results are shifted a little bit upwards compared to 
the exact results, but this error is an order of magnitude 
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FIG. 15: Exact KS (middle) and TDDFT-corrected (ALDA) 
results as a function of the wall radius. The singlet is above, 
and the triplet below, the KS curve.. 



smaller than those due to approximate kernels. 

We conclude that using EXX potentials is clearly suf- 
ficient for these calculations. On the other hand, the 
LDA results diverge. Thus, without the special handling 
of Ref. (42|, standard approximations (LDA, GGAs, hy- 
brids) arc useless here. 

D. Effect of the wall 

In this section, we demonstrate our claim that the in- 
finite wall can be placed quite close to the atom, as long 
as the part of the potential that differs from — 1/r is es- 
sentially zero beyond the wall (Eq. (|2.4|) ). In Fig. [151 wc 
show the singlet and triplet s-wave phase shifts for wall 
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FIG. 16: Same as Fig. 1121 for p-wave. Color figure online. 



FIG. 17: Same as Fig. 1131 for p-wave. Color figure online 



distances between 5 and 8 bohr, and compare the results 
with i?b = 100 a.u. The smallest distance is determined 
by the requirement that the box does not significantly 
distort the self-consistent potential in the region of signif- 
icant density, i.e., occupied orbitals must be unaffected. 
The results for the small wall distances indeed correspond 
to the correct phase shifts. By varying the wall distance 
at small Rb we can fill out the complete curve, giving an 
alternative to performing a single calculation with a very 
large Rb- This is an important conclusion if one wants 
to use finite basis-set methods instead of a grid based 
code [48jj, because the range of finite basis-sets, even if 
they include very diffuse functions, is limited to much 
shorter radii. 
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FIG. 18: s-QD and s-wave phase shift for Hartree and ALDA 
compared to accurate data [Hoi ]. The bottom three curves 
correspond to singlet, the top three to triplet. 



E. Different /-values 



Until now we focussed on s-wave scattering and s-QDs. 
Of course we can also study higher /-values, using the /- 
dependent Eq. (|2.12p for the phase shift. We show the 
difference between our p-QD and phase shift (i.e. / = 1) 
and the reference values for He in Figs. [16] and [171 In 
Rcf. [49j . we also showed results for e?-wavc scattering. 

Adding xc-kernels notably improves the Hartree results 
just as in case of the s-QD and s-wave phase shift. For 
negative and low energies, the hybrid performs the best of 
the three kernels, but for larger energies ALDA performs 
better, being very close to the reference results. The bare 
EXX kernel does not perform as well in this case, unlike 
for s-wave phase shifts. Thus ALDA is better overall, 
and can easily be applied to any system tractable with 
ground-state DFT. 



F. Scattering from a neutral system: e-H scattering 

When we scatter from a positive ion like He + as we 
did throughout this paper, the (JVt + l)-electron system 
is a neutral atom with a long-ranged KS potential. In 
Ref. [IH we studied another two-electron system, namely 
electron scattering from the hydrogen atom. In this case 
the (Nt + l)-electron system, H~, is a negative ion with 
a short-ranged KS potential. Because of this, there is no 
Rydberg series or QD and only one bound state. We now 
calculate the phase shift using Eq. (|2.10|) . i.e. relative to 
free potentials. 

We show our results in Fig.[l8j We see that, just as for 
e-He + scattering, the KS result, which is equal to triplet 
Hartree, lies in between both the accurate singlet and 
triplet values. ALDA corrects this and the ALDA results 
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are again very close to the reference data. So TDDFT 
works well for electron scattering from both cations and 
neutral atoms. For more details on this particular system 
we refer to Ref. fl9l ]. 

VI. CONCLUSIONS 

Electron scattering from H (as shown in Fig. [T5|) and 
electron scattering from He + (as outlined in this paper) 
are two-electron systems. But our TDDFT methods can 
also be applied to elastic electron-scattering for systems 
with more than two electrons. 

Results for electron scattering from Li can be found 
in Ref. (49[. The (Nt + l)-electron system is Li - and 
has a short-ranged potential. The singlet TDDFT p- 
wave phase shift in this system suffers from an artifact 
that appears very close to zero energy (energies smaller 
than 1 eV). But for higher energies the TDDFT values 
correspond very well to the available reference data, as 
do the triplet p-wave and s-wave phase shifts. 

A many-electron system with a long-ranged potential is 
Be, corresponding to electron scattering from Be + . Scat- 
tering phase shifts of e-Be + can be found in Ref. [Blj . 
In this system there are a number of low-energy reso- 
nances of double excitation character, not captured by 
our TDDFT method. Such effects require going beyond 
the adiabatic approximation of TDDFT and using a fre- 
quency dependent kernel [52l |. But apart from missing 
these resonances TDDFT again performs well (HT| ■ 

So in general this method is expected to work well for 
low-energy electron-atom scattering problems, without 
resonances due to multiple excitations, where the (Nt + 
l)-electron system is bound. We have not yet addressed 
the complications of non-spherical systems or of open- 



shells. 

Another point we have not addressed are the oscillator 
strengths that correspond to the TDDFT excitation en- 
ergies and how they behave when moving across the con- 
tinuum. This situation is studied in Ref. [53j |. It is found 
that the KS oscillator strength corresponding to the ion- 
ization energy (exactly given by minus the energy of the 
KS HOMO) does not correspond to the exact oscillator 
strength. But the exact and KS oscillator strengths do 
merge at higher energies and several results for 2-electron 
systems are given in the reference. 

In this paper we reviewed our theory and calculations 
for obtaining both the quantum defect and electron-atom 
scattering phase shifts using TDDFT. We have shown 
that our method of putting the system in a finite hard 
wall cavity in order to extract the continuum states and 
phase shifts works within TDDFT. In our previous work 
we showed that the method works well for the short- 
ranged potential system e-H _ . In this paper we first 
illustrated our box method on a (relatively) simple toy 
problem, in which the answer was known exactly. Next, 
we showed that the method works equally well for the 
long-ranged helium atom and e-He + scattering. We also 
showed that the TDDFT quantum defect and phase shifts 
transition smoothly across the zero energy boundary, in 
accordance with Seaton's theorem. 

Our final conclusion is that TDDFT works well for 
high- lying excitations in neutral atoms and cations. To 
infinity, and beyond! 
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